000D0 http:Wwww.cqvip.com 


动物 学 研究 2005, Jun. 261(3): 225 -229 CN 53-1040/Q ISSN 0254 - 5853 
Zoological Research 


一 种 在 核 并 酸 水 平 检测 自然 选择 的 新 方法 


2 
(1. 西北 农林 科技 大 学 动物 科技 学 院 , 陕西 杨陵 “712100; 2. 中 国 科学 院 昆明 动物 研究 所 细胞 与 分 子 进化 重点 实验 室 ， 
中 德 马 普 青 年 科学 家 小 组 ,云南 昆明 650223; 3. 徐州 师范 大 学 生物 技术 研究 所 ， 江 苏 徐州 ”221116) 


摘要 : 非 编码 区 序列 在 基因 表达 调控 中 起 着 重要 作用 ,但 其 在 进化 过 程 中 是 否 受 到 选择 作用 一 直 较 难 检 
测 。 最 近 有 一 些 研究 使 用 平均 的 核 苷 酸 蔡 换 速率 与 中 性 序列 的 核 苷 酸 蔡 换 速率 的 比值 (。) 作为 检测 非 编码 区 
总 体 受 选择 作用 的 指标 ; 但 是 对 于 非 编码 区 而 言 ， 了 解 具 体 哪些 核 背 酸 受 到 选择 作用 更 具有 意义 。 我 们 借鉴 
Nielsen & Yang (1998) 检测 单个 氨基 酸 位 点 是 否 受 选择 作用 的 思路 ， 在 最 大 似 然 法 的 模型 下 ， 提 出 一 种 在 核 苷 
酸 位 点 水 平 上 对 自然 选择 作用 检测 的 方法 。 本 方法 能 够 检测 在 进化 过 程 中 对 功能 分 化 有 重要 贡献 的 核 苷 酸 位 
点 ， 包 括 编码 和 非 编码 区 。 将 此 方法 应 用 于 熟知 的 受到 正 选 择 作 用 的 蛋白 编码 基因 序列 ( HIV-1 包装 蛋白 基因 
编码 区 ) ， 均 能 够 检测 到 那些 已 知 的 受到 正 选 择 的 核 背 酸 (密码 子 ) 位 点 ， 说 明 此 方法 可 以 有 效 地 在 核 苷 酸 位 
点 水 平 检测 选择 作用 ; 又 将 此 方法 应 用 于 非 编码 区 (CTCF 基因 5'UTR)， 也 得 到 了 良好 的 结果 。 
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Abstract: Although noncoding regions play an important role in gene expression and regulation, it is difficult to de- 
tect natural selection at this level. Recently, some studies use the ratio (w) between the nucleotide substitution rate in the 
detected regions and the nucleotide substitution rate in neutral regions as an indicator to detect natural selection in noncod- 
ing regions. However, to the noncoding regions ，it's more informative to identify those nucleotide sites under positive se- 
lection. We developed a new maximum-likelihood method to detect natural selection at the nucleotide site level and 
to identify those nucleotide sites that may contribute to functional divergence. This method can be applied to both coding 
regions and noncoding regions. Appling this method to previous reported genes that subjected to positive selection shows 
that this method is efficient to detect natural selection on nucleotide sites both in coding regions and noncoding regions . 
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在 分 子 水 平 检 测 选 择 作用 是 近年 来 进化 遗传 学 。 有 蛋白 的 基因 ， 错 义 替换 率 (nonsynonymous substitu- 
中 的 一 个 热门 话题 (Tajima, 1989; Fu & Li, 1993; tion rate，Ka) 与 同 义 替 换 率 (synonymous substitu- 
Hudson et al, 1987; Mcdonald & Kreitman, 1991; tion rate，Ks) 的 比值 (Ka/Ks) 被 广泛 地 用 来 检测 
Fay & Wu, 2001; Zhou & Wang，2004)。 对 于 编码 ” 基因 是 否 受 到 选择 作用 (Miyata & Yasunaga， 
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1980)。 在 Ks 为 中 性 替换 的 假设 下 ， 如 果 Ka/Ks < 
1, 说 明 此 基因 在 进化 过 程 中 受到 纯化 选择 (puri- 
fying selection) 的 作用 ; 若 Ka/Ks = 1， 即 同 义 替换 
与 错 义 替换 有 相同 的 可 能 被 固定 下 来 ， 指 示 此 基因 
经 历 的 是 一 个 中 性 的 进化 过 程 ， 而 如 果 Ka/Ks > 1， 
即 某 些 错 义 替换 具有 选择 上 的 优势 ， 其 固定 的 概率 
和 速度 要 大 于 同 义 替换 ,说明 此 基因 在 进化 过 程 中 
受到 了 正 选 择 (positive selection )。 我 们 在 证 明 与 
雄性 生殖 相关 的 基因 受到 强烈 正 选 择 作用 时 即 采 用 
了 这 一 思路 (Wyckoff et al，2000)。 然 而 由 于 在 进 
化 过 程 中 ， 一 个 基因 的 大 部 分 位 点 在 大 部 分 时 间 
中 ， 一 般 被 认为 是 受到 强 的 功能 限制 ， 因 此 用 整体 
基因 的 平均 Ka/Ks > 1 作为 受到 正 选 择 的 标准 由 于 
过 于 严格 而 缺乏 检测 效力 (Yang & Bielawski， 
2000; Yang，2002b; Endo et al ,1996)。 随 后 人 们 采 
用 分 割 法 ,将 整个 基因 按 功 能 分 成 几 个 区 域 分 别 计 
算 各 个 区 域 的 Ka/Ks， 从 而 提高 检验 的 效力 
(Hughes & Nei，1988)。 然 而 与 功能 进化 密切 相关 
的 位 点 在 DNA 序列 上 并 不 一 定 是 连续 排列 的 ， 因 
此 这 种 改进 对 于 检验 效力 的 提高 仍然 是 有 限 的 。 最 
近 有 研究 者 提出 一 系列 基于 似 然 法 和 简约 法 针对 单 
个 氨基 酸 位 点 进行 选择 作用 检测 的 统计 学 方法 ， 这 
些 方法 不 但 能 够 更 有 效 地 检测 基因 在 进化 历程 中 是 
否 受 到 选择 的 作用 ， 并 且 可 以 预测 出 那些 在 进化 过 
程 中 对 功能 分 化 有 重要 贡献 的 、 受 正 选 择 作 用 的 氨 
基 酸 位 点 (Nielsen & Yang, 1998; Fitch et al, 1997; 
Suzuki & Gojobori, 1999; Yang & Nielsen ，2002 )。 
因此 这 些 方法 在 很 大 程度 上 提高 了 Ka/Ks 检验 的 检 
测 效力 ， 使 得 对 基因 适应 性 进化 的 检测 由 整体 水 平 
或 是 区 域 水 平 还 原 到 单个 氨基 酸 位 点 的 水 平 上 。 
相对 于 编码 蛋白 基因 分 子 进化 的 大 量 研究 ， 人 
们 对 非 编 码 区 ， 包括 RNA 基因 、UTRs (untranslat- 
ed regions) 、 启 动 子 、 内 含 子 和 基因 间 序 列 的 进化 
研究 远 远 不 够 。 近 年 来 越 来 越 多 的 研究 也 表明 非 编 
码 区 同样 存在 着 选择 作用 (Makalowski et al，1996; 
Makalowski & Boguski, 1998; Larizza et al, 2002; 
Michael et al，2004)。 在 对 非 编码 区 的 分 析 中 ， 类 
似 于 前 面 提 到 的 Ka/Ks 的 标准 ， 被 检测 区 域 核 苷 酸 
的 替换 速率 (Kd) 与 中 性 区 域 的 核 苷 酸 替 换 速率 
(Kn) 的 比值 。(w = Kd/Kn) 也 被 用 来 作为 这 些 
区 域 所 受 选 择 作 用 的 指标 。w <1, w=1 和 ww>1 
分 别 作 为 纯化 选择 ， 中 性 进化 和 正 选 择 的 指示 。 在 
我 们 研究 一 个 年 轻 的 非 蛋 白 编码 的 RNA 基因 





sphinx 基因 时 ， 我 们 使 用 类 似 的 思路 证 明了 sphinx 
基因 受到 了 正 选 择 的 作用 (Wang et al, 2002)。 但 
是 , 与 前 面 检 测 编码 蛋白 基因 所 存在 的 问题 相似 ， 
这 种 整体 替换 速率 的 比较 仍然 存在 检验 效力 较 低 的 
问题 ， 通 常 只 有 在 选择 作用 很 强 时 才能 检测 到 。 为 
了 解决 这 个 问题 ， 本 文 尝试 将 Nielsen & Yang 
(1998) 的 方法 从 氨基 酸 位 点 的 水 平 上 推广 到 单个 
核 苷 酸 位 点 的 水 平 上 ， 这 样 ， 此 方法 不 仅 能 够 检验 
编码 蛋白 基因 的 选择 作用 ， 更 重要 的 是 它 也 能 够 对 
非 编码 区 的 选择 作用 进行 检测 。 


1 模型 与 方法 


1.1 位 点 的 分 类 
考虑 多 条 已 排序 的 核 苷 酸 序 列 ， 并 且 这 些 序列 
间 的 系统 发 育 关系 已 知 。 假 设 序列 上 各 个 位 点 的 进 
化 相互 独立 。 我 们 将 所 有 的 核 苷 酸 位 点 分 为 3 类 : 
中 性 位 点 、 保 守 位 点 和 正 选 择 位 点 。ri、rz 和 rs 
分 别 表示 待 检测 序列 的 中 性 位 点 、 保 守 位 点 和 正 选 
择 位 点 的 替换 速率 ，ro 表示 中 性 序列 的 替换 速率 ， 
它们 的 替换 速率 与 中 性 序列 的 替换 速率 的 比值 定义 
如 下 : 
wi! = riro = 1 
w2 = r2/ro = 0 
w3 = r3/ro > 1 
我 们 将 保守 位 点 和 中 性 位 点 的 替换 速率 定 为 固 
定 值 0 或 1 是 为 了 计算 方便 ， 也 可 以 不 限制 wm! 和 
w2 的 值 而 将 其 作为 未 知 变量 用 似 然 法 估计 而 使 模 
型 更 符合 实际 。 
1.2 中 性 替换 速率 的 计算 
考虑 图 1 所 示 由 4 条 序列 构成 的 系统 发 育 树 。 
我 们 根据 与 这 4 条 待 检 验 序列 相关 的 中 性 序列 如 内 
含 子 、 假 基因 、 基 因 间 序列 ， 或 者 同 义 替换 位 点 ， 
可 用 PAML 或 PHYLIP 程序 中 的 最 大 似 然 法 估计 各 
个 枝 长 (Yang, 2002a; Felsenstein, 2002)。 这 样 得 
到 6 条 枝 的 枝 长 ， 所 得 枝 长 即 为 各 个 枝 的 中 性 替换 
速率 ， 分 别 表示 为 Vis，V2s，.. 
1.3 最 大 似 然 法 估计 未 知 参 数 
为 方便 计算 ， 我 们 仅 考 虑 Felsenstein 的 核 苷 酸 
替换 模型 (Felsenstein，1981 ) ， 此 方法 不 考虑 转换 
颠 换 的 差异 ， 而 根据 序列 中 某 种 核酸 的 比例 确定 向 
其 突变 的 概率 。 其 他 替换 模型 计算 方法 类 似 ， 只 需 
添加 更 多 的 待 估 参 数 即 可 。 由 于 我 们 事先 并 不 知道 
每 个 位 点 属于 哪 一 类 位 点 ,因此 假设 整 条 序列 中 ， 
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图 1 四 条 序列 构成 的 系统 发 育 树 
Fig. 1 Phylogenetic tree of 4 sequences 


中 性 位 点 、 保 守 位 点 和 正 选 择 位 点 所 占 的 比例 分 别 
为 pl，p2，p3，p1+ pz+p3=1。 我 们 采用 最 大 似 
然 法 估计 3 个 未 知 参数 (wa, pl, Pp2)o 假设 这 4 条 
序列 长 度 为 n， 考 虑 其 中 的 位 点 x;， 根 据 现在 观察 
到 的 4 条 序列 在 此 位 点 上 的 数据 ， 它 的 似 然 值 为 


f(xi) = Dn x f(xi | rm) ，pi 为 此 位 点 属于 第 


类 位 点 的 概率 ， 了 (xilrs) 表示 当 此 位 点 属于 第 上 
类 位 点 时 ， 出 现 现 有 情况 的 条 件 概率 ， 计 算 它 时 要 
考虑 每 一 条 枝 的 进化 情况 。 对 于 图 1 所 示 情 况 来 说 
(xilre) = PPzP5Pe Psx6Pk6， 户 表示 某 位 点 
上 由 结 点 i 上 的 核 苷 酸 蔡 换 为 结 点 上 的 核 苷 酸 的 
概率 。 根 据 Felsenstein (1981) 的 核 车 酸 替 换 模型 : 
Ti+ (ln) x exp(-— wiVy/(1 — 4 )), 当 结 
p， -| 点 上 与 结 点 i 上 的 核 车 酸 相 同时 ; 
” mm-mxexp(- wrVy/(1 -5)), 当 结 点 i 与 
结 点 六 上 的 核 苷 酸 不 同时 。 
0 


j=C，T，C 或 者 A; b=》 。 这 样 便 计算 出 了 
第 i 个 位 点 的 似 然 值 。 同 理 ， 对 于 个 位 点 都 做 此 计 


算 ， 求 它们 的 对 数 和 ， 记 为 工 , 了 = Pin) , 


L 即 为 整 条 序列 总 的 似 然 值 。 调 整 w3, pi, Pp2 各 参 
数值 使 荆 最大， 估计 出 w3, pi, pz。 
1.4 似 然 比 检验 

以 上 是 考虑 序列 中 含有 正 选择 位 点 的 情况 。 同 
理 我 们 也 考虑 假设 序列 中 不 含 正 选择 位 点 的 情况 ， 
即 序列 仅 由 中 性 位 点 和 保守 位 点 组 成 ， 仍 然 可 以 计 
算出 一 个 似 然 值 ， 记 为 L1， 这 是 在 中 性 模型 下 计 
算 的 似 然 值 。 然 后 比较 工 与 是 否 有 显著 差别 ， 


即 检 验 用 选择 模型 是 否 能 比 中 性 模型 更 好 地 解释 数 
据 。 由 于 选择 模型 包含 了 中 性 模型 ， 且 比 其 多 2 个 
未 知 参数 ， 我 们 采用 似 然 比 检验 ， 即 2 ( 工 - 工 ) 
近似 服从 自由 度 为 2 的 x? 分 布 。 
1.5 后 验 概率 的 计算 

当 上 述 的 检验 显著 后 ， 我 们 可 以 计算 每 个 位 点 
属于 某 一 类 位 点 的 后 验 概率 ， 即 Pyos = ps x (xil 
rk)，Ppost 表 示 第 i 个 位 点 属于 第 类 位 点 的 后 验 概 
率 。 若 某 位 点 属于 第 三 类 位 点 的 概率 很 大 ， 假 如 大 
于 0.95 或 0.99， 我 们 就 认为 此 位 点 在 进化 过 程 中 
受到 正 选 择 。 


2 结果 与 分 析 


2.1 对 蛋白 编码 基因 的 分 析 

基于 以 上 方法 ,我 们 首先 对 HIV-1 的 包装 蛋白 
基因 做 分 析 ， 由 于 没有 内 含 子 序 列 ， 我 们 用 同 义 替 
换 位 点 作为 中 性 序列 来 估计 中 性 替换 速率 。HIV-1 
包装 蛋白 基因 序列 数据 取 自 Nielsen & Yang 
(1998) ， 计 算 结果 见 表 1。 


表 1 HIV-1 包装 蛋白 基因 分 析 结果 
Tab. 1 Testing result on the HIV-1 envelope gene 


pl Pp2 Ww3 L Li 


0.2 0.48 10 -546.762 2 -3567.568 1 


用 工 和 万 进行 似 然 比 检验 ， 即 检验 是 否 选择 
模型 能 更 好 地 符合 已 知 数据 ， 检 验 结果 极 显 著 〈( 己 
<0.01)， 说 明 此 基因 中 存在 一 定 比 例 受 到 正 选 择 
的 位 点 。 

对 每 个 位 点 计算 其 属于 正 选择 位 点 的 后 验 概 
率 ， 结 果 显示 有 4 个 位 点 的 后 验 概 率 值 极 显著 〈 表 
2)。 同 时 用 Nielsen & Yang (1998) 的 方法 对 序列 
进行 分 析 ， 亦 检测 出 4 个 极 显著 位 点 ， 其 中 3 个 位 
点 完全 相同 ，1 个 位 置 上 比较 接近 。 说 明 我 们 的 方 
法 能 有 效 地 检测 受到 正 选 择 的 核 车 酸 位 点 ( 表 2)。 
2.2 对 非 编码 序列 的 分 析 

由 于 本 方法 能 够 直接 检测 核 苷 酸 位 点 ， 因 而 可 
以 对 非 编 码 区 进行 分 析 。 我 们 以 CTGF 基因 (con- 
nective tissue growth factor) 为 例 。 此 前 ，Larizza et 
al (2002) 用 人 、 鼠 和 牛 中 的 CTGF 基因 的 mRNA 
序列 ， 基 于 对 整个 序列 的 核 苷 酸 替 换 率 与 同 义 替 换 
率 的 比较 ， 认 为 此 基因 的 5'UTR 在 进化 过 程 中 受到 
正 选 择 的 作用 ,我 们 用 本 方法 也 对 此 基因 的 5'UTR 
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表 2 HIV-1 包装 蛋白 中 预测 的 正 选择 位 点 
Tab. 2 Predicted positively selected sites in the HIV-1 envelope gene 
本 文 方法 Method of this study Yang 的 方法 Method of Yang* 

位 点 所 在 的 核 苷 酸 位 置 位 点 对 应 的 氨基 酸 的 位 置 后 验 概率 位 点 所 在 的 氨基 酸 位 置 后 验 概率 
Position of predicted sites Position of predicted sites we Position of predicted sites Sa 
， 5 : 了 Posterior probability Posterior probability 
in the nucleotide sequence in the protein sequence in the protein sequence 

83 28 0.991 624 1 28 0.992 6 

107 36 0.991 472 2 31 0.982 1 

198 66 0.991 1107 66 0.992 5 

203 68 0.991 400 7 68 0.995 5 


”See Nielsen & Yang (1998). 
表 3 CTGF 基因 分 析 结 果 
Tab. 3 Testing results on the 5’UTR of the CTGF gene 


pi Pp2 WwW3 L Li 
0.1226 0.3724 70.5 -251.0027 -271.6435 





表 4 CTGF 基因 5'UTR 中 预测 的 受到 正 选 择 作用 


的 位 点 
Tab. 4 Predicted positively selected sites in the 5’ UTR 
of CTGF 
位 点 所 在 的 核 苷 酸 位 置 ( 距 翻 
译 起 始 位 点 上 游 的 碱 基 数 ) 
Se of Predicted sites in 后 验 概率 Posterior probability 
e nucleotide sequence 
(Nucleotide numbers from 
upstream of start coden ) 
4 0.973 694 9 
12 0.965 945 4 
14 0.983 540 1 
34 0.955 975 7 
38 0.951 803 4 
40 0.991 8306 
55 0.974 840 9 
73 0.983 540 1 
74 0.976 646 7 
77 0.956 469 1 
90 0.968 406 
112 0.965 945 4 
114 0.968 406 
115 0.975 347 3 
117 0.976 646 7 


做 了 分 析 ， 各 参数 的 估计 值 见 表 3， 选 择 模 型 与 中 
性 模型 的 似 然 比 检验 极 显著 (P <0.01), 说 明 此 基 
因 的 5'UTR 区 的 确 存 在 正 选择 作用 ， 并 检测 出 一 些 
后 验 概率 显著 的 位 点 〈 表 4)。 更 多 的 近 缘 物 种 间 该 
基因 序列 的 比较 或 人 工 突变 研究 可 以 进一步 证 实 这 
些 结论 。 


3 讨论 
3.1 对 蛋白 编码 基因 选择 作用 分 析 的 有 效 性 


由 对 HV-1 包装 蛋白 基因 的 分 析 可 以 看 出 ， 
本 方法 所 得 结果 与 以 氨基 酸 为 分 析 单 位 的 Nielsen & 
Yang (1998) 的 方法 所 得 结果 基本 一 致 ， 说 明 对 于 
编码 蛋白 基因 ， 本 方法 是 有 效 的 。 事 实 上 ， 如 果 以 
同 义 替 换 速 率 作 为 中 性 替换 速率 ， 那 么 在 错 义 替换 
率 大 于 同 义 替 换 率 的 氨基 酸 位 点 (密码 子 ) ， 其 平 
均 核 苷 酸 蔡 换 速 率 肯 定 大 于 中 性 替换 速率 。 因 此 ， 
这 两 种 方法 在 本 质 上 是 相似 的 ， 应 当 得 到 基本 一 致 
的 结果 。 

然而 ， 由 于 密码 子 使 用 偏好 等 原因 ， 同 义 替 换 
速率 未 必 是 中 性 的 (Akashi，1994)。 因 此 ，Ka/Ks 
大 于 1 既 可 能 是 由 于 正 选 择 所 驱动 的 快速 的 错 义 替 
换 率 ， 也 可 能 是 由 于 同 义 替换 位 点 受到 选择 限制 所 
导致 。 所 以 用 Ka/Ks 是 否 大 于 1 作为 是 否 受 到 正 选 
择 的 标准 可 能 会 引起 错误 的 判断 。 我 们 在 分 析 中 采 
用 内 含 子 、 假 基因 或 基因 间 序 列 作为 中 性 序列 时 可 
以 尽量 避免 这 种 错 判 。 尽 管 寻找 一 条 绝对 的 中 性 序 
列 是 很 困难 甚至 是 不 可 能 的 ， 即 便 是 内 含 子 等 人 们 
一 般 认 为 的 中 性 序列 也 可 能 受到 选择 作用 (Parsch， 
2003) ,但 是 各 种 类 型 的 序列 受到 选择 作用 影响 的 
程度 是 不 同 的 ， 因 此 尽 可 能 准确 地 估计 中 性 替换 速 
率 是 本 方法 得 到 较 可 信 结 果 的 关键 。 我 们 认为 ， 使 
用 待 检 序列 邻近 的 中 性 序列 应 是 较 好 的 方法 。 
3.2 对 非 编码 区 选择 作用 分 析 的 有 效 性 

相对 于 Nielsen & Yang (1998) 的 方法 ， 本 方法 
的 最 大 优点 就 是 可 以 对 非 编 码 区 进行 分 析 ， 对 CT- 
GF 基因 5'UTR 的 分 析 结 果 说 明了 其 有 效 性 。 由 于 
非 编码 区 被 认为 对 基因 的 表达 调控 有 重要 作用 ， 因 
此 该 区 域 所 受 的 选择 作用 有 时 对 基因 的 进化 是 至 关 
重要 的 。 然 而 如 何 确定 非 编 码 区 中 的 功能 序列 历来 
是 研究 的 一 个 难题 ， 通 过 本 方法 预测 的 受到 正 选 择 
的 位 点 可 以 为 进一步 的 功能 研究 提供 一 个 理论 参 
考 。 
3.3 存在 的 问题 
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由 于 核 苷 酸 的 种 类 只 有 4 种 ， 因 此 对 于 n 条 已 
对 齐 (alignment) 的 序列 来 说 ， 每 个 核 苷 酸 位 点 理 
论 上 可 以 有 4" 种 可 能 。 而 组 成 蛋白 质 的 氨基 酸 有 
20 种 ， 编 码 20 种 氨基 酸 的 密码 子 有 61 个 ， 因 此 ， 
以 氨基 酸 (密码 子 ) 为 分 析 单 位 时 ， 会 利用 3 个 核 
苷 酸 位 点 的 信息 ， 其 每 个 氨基 酸 位 点 理论 上 有 61" 
种 可 能 。 这 样 对 于 编码 蛋白 基因 来 说 ， 这 种 基于 核 
苷 酸 水 平 的 分 析 方 法 与 以 氨基 酸 为 分 析 单 位 的 方法 
相 比 ， 对 于 区 分 不 同位 点 间 的 选择 作用 分 辨 率 相对 
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